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ABSTRACT 

Context. With the development of new polarimeters for large telescopes, the spectro-polarimetric study of astrophysical bodies is 
becoming feasible and, indeed, more frequent. In particular, this is permitting the observational study of stellar magnetic fields 
Aims. With the aim to optimize and interpret this kind of observations, we have produced a spectral atlas of circular polarization 
in a grid of stellar atmospheric models with effective temperatures between 3500 and 10000 K, surface gravities log(g) = 3.5 - 5, 
metallicities between 10~ 2 and 1, and magnetic field strengths of 100, 1000 and 5000 G 
Methods. We have computed the emergent Stokes / and V flux spectra in LTE of more than 10 5 spectral lines 

Results. The atlas and several numerical tools are available in electronic format and may be downloaded from 
http://download.hao.ucar.edu/pub/PSA/ In this paper we review and discuss some of its most relevant features, such as which spectral 
regions and individual lines harbor the strongest signals, what are interesting lines to observe, how to disentangle field strength from 
filling factor, etc. 
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1. Introduction 

The development of large-aperture telescopes and high- 
sensitivity polarimetric instrumentation is providing the 
Astrophysical community with new tools to explore the uni- 
verse. Polarimetric observations are opening new and exciting 
possibilities for the diagnostics of light-matter interaction 
processes in many different scenarios (see, e.g., the various 
reviews in ? ?). 

There are two major reasons why spectral lines produced in 
stellar atmospheres may become polarized. First, we may have 
scattering processes in which the photons emitted by the source 
are deflected towards Earth. Depending on the details of the in- 
teraction, one receives light that is (at least in part) linearly po- 
larized in a direction that is closely related to the scattering ge- 
ometry (and some times other factor, such as the magnetic field). 
This is obviously of great interest for the detection and charac- 
terization of exoplanets (e.g., ? ?), protoplanetary and accretion 
disks (e.g., ??;??;? ?), binary systems (e.g., ? ?) etc. 

Another source of polarized light is the Zeeman effect that 
splits the atomic energy sublevels when magnetic fields are 
present in the spectral line formation region. The strength of 
the magnetic field required to produce a certain amount of po- 
larization depends on details of the radiative transfer (Zeeman 
sensitivity of the spectral lines, opacity profile, source function 
gradient, etc). Angular integration of a spherically-symmetric 
emitting source results in circular polarization spectral profiles 
(linear polarization would cancel out). Analysis of these spectral 
profiles provides unvaluable clues to study, e.g. magnetic activ- 
ity in stars (e.g., ? ?; ? ?). 

For several decades now, the grid of radiative equilibrium 
model atmospheres and spectra computed by Kurucz (for an up- 



dated reference see ? ?) have been very useful in the preparation 
and interpretation of stellar observations. We felt that a polariza- 
tion atlas of the signal emergent from those models, for several 
magnetic field strengths and metallicities, would also be useful 
for future observations of stellar magnetic fields. The entire grid 
of polarized spectra computed in this work is publicly available 
for download from the following URL: 

\protect\vrule widthQpt\protect\href {http : //download . hao . u 



2. The Atlas 

The synthetic spectra have been obtained using a compilation of 
specttal line data from ? (?) between 400 and 1600 nm for the 
model with r e ff=5750 K and log(g)=4.5 (hereafter, the solar- 
like model), and between 400 and 1000 nm for the rest (if there 
is enough interest, we may consider extending the wavelength 
range). This is a useful range to consider, as it is accessible to 
modern spectro-polarimettic instruments such as MUSICOS (? 
?) or ESPaDOnS (? ?). 

Due to insufficient data on partition functions we restricted 
ourselves to lines of up to the third ionization stage. Lines of 
higher ionization stage have been excluded. The effective Lande 
factor (g e ff) for lines that are not well described by the LS cou- 
pling or for which there is no level configuration information 
has been set to zero, so no polarization signal will appear in 
these linefl Following ? (?), the Ha line is assumed to have 
an effective Lande factor of 1 . A comprehensive list of the lines 
excluded or computed with zero Lande factor (35042 for the 
solar-like model and 27188 for the others) is available for down- 
load in electronic format. The total number of lines considered 
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1 Although we shall often refer to effective Lande factors, our com- 
putation actually takes into account the full Zeeman pattern (anomalous 
splitting) and does not resort on the triplet approximation. 
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is 165501 in the extended range for the solar-like model and 
128794 for all others. Table Q] lists the distribution of spectral 
lines included and excluded in each wavelength range. 

The full polarized radiative transfer vector equation for a 
given wavelength is usually written as: 



Table 1. Number of lines included and excluded in the syntheses 



dl(s) 

ds 



-K(s)[I(s)-S(s)]. 



(1) 



where s is the distance along the ray path, I is the Stokes vec- 
tor, S is the source function vector (which in LTE is usually the 
Planck function times the transpose of the unitary [1,0,0,0] vec- 
tor), and finally K is the absorption matrix: 



K 



( m t]Q rju rj v 

t]q m pv -rjv 

vu -pv m p Q 

V rj V pu -p Q i], 



(2) 



For details on the meaning of the symbols employed in the equa- 
tions above, the reader is referred to ? (?). Note that, although we 
are only interested on the emergent intensity / and circular po- 
larization V, we must compute all the four Stokes parameters 
consistently due to the off-diagonal terms in K. 

The syntheses were carried out using a modified version of 
the code LILIA (? ?). This code was designed for the synthe- 
sis and inversion of polarized atomic spectral lines formed in 
LTE (molecular bands cannot be treated with this code). Solar 
metal abundances are those from ? (?). For all other metallic- 
ities, abundances are scaled by a constant factor. Background 
opacities are obtained with the opacity package used by ? (?). 
The atmosphere is assumed to be one-dimensional and plane- 
parallel. The formal solution of the polarized transfer equation is 
based on the Hermitian method of ? (?). However, a much sim- 
pler method may sometimes be used (see the Weakly Polarizing 
Media approximation below). 

The plane-parallel approximation is justified by our choice 
of not considering the largest stellar models (log(g) between 3.5 
and 5). The LTE approximation implicit in the expressions of 
the source function and the absorption matrix is justified by the 
range of effective temperatures, T e s between 3500 and 10000 K. 

For consistency, we started with the mass scale from the orig- 
inal models and transformed it into an optical depth scale, using 
our opacities in the calculation. The optical depths obtained are 
the reference height scale for the radiative transfer. In order to 
avoid very hot (where LTE may not be a good approximation for 
many photospheric lines) or very large stars (where the spheri- 
cal geometry may be needed), we decided to restrict ourselves to 
a range of effective temperatures (hereafter T e g) between 3500 
and 10000 K and surface gravities (hereafter log(g)) between 
10 3 5 and 10 5 . The mixing-length treatment of these models has 
been described by ? (?) (see also ? ?). 

We selected five different representative metallicities, with 
abundance ratios (relative to solar) of 1, 10~ 5 , 10 -1 , 10~' 5 , and 
10~ 2 . For each abundance, T e s and log(g), we computed three 
different magnetic field strengths of 100 G, 1000 G and 5000 G. 

It is clear that the amount of all possible realizations of mag- 
netic field configurations on the surface of a star cannot be ex- 
hausted by any reasonable grid of models; that is certainly not 
the aim of this work. To find out the particular configuration that 
best fits some observation, one would need to resort on least- 
squares inversion methods (such as that of ? ?). Our atlas has 
been calculated to give a qualitative account of actual observa- 
tions, as well as to explore possible observational strategies, find 
the most promissing spectral lines for magnetic field diagnostics, 
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etc. Therefore, we have assumed the simplest possible configura- 
tion, given by a monopolar magnetic field pointing radially away 
from the center and of constant strength. 

The spectra are calibrated in flux using a 3-point Gaussian 
quadrature. We looked at several randomly-chosen models and 
verified that increasing the quadrature grid to 4 or 5 points did 
not make any significant difference. 

The models used here have a depth-independent microturbu- 
lent velocity of 2 km s _1 and a macroturbulence of 1 km s _1 . 
These are roughly representative solar values . In addition to tur- 
bulence, other line-broadening mechanisms are taken into ac- 
count after ? (?): 

- Radiative broadening. We set F,. = 0.22233 /A 2 (with A in 
cm). This formula is only a rough approximation but radia- 
tive broadening is usually very small anyway. 

- Van der Waals broadening. Note that the ? (?) expressions 
are only accurate for neutral or singly ionized atoms. For 
other ions we neglect this broadening. 

- Stark broadening. We use the approximation Fs -10.06666 
+ log 10 (P c )-5/6 log 10 (r), where P e and T are the electron 
pressure and the local temperature (in c.g.s. units). 

Even assuming LTE, the synthesis of the Stokes profiles of 
over 10 5 spectral lines in a large grid of model atmospheres is 
a very demanding computing task. To bring the computing time 
within reasonable practical limits we implemented a number of 
optimizations that do not harm significantly accuracty of the po- 
larized spectra: 

- Background opacities. These are not recomputed at each 
wavelength. Instead, they are evaluated every 0.1 nm. This 
is noticeable, e.g. when one inspects closely the continuum 
intensity over a very short wavelength range. FigureQ]shows 
the continuum for the solar-like model. A very close zoom- 
in reveals the discretization in the continuum caused by this 
approximation. Notice, however, the ordinates scale which 
shows that the continuum is virtually flat and thus the error 
introduced (the "step" of the staircase) is very small, at the 
10~ 4 level. 

- Weak lines. Depending on the particular physical conditions 
on a given atmosphere, many spectral lines turn out to be 
negligibly weak. Removing these lines from the computation 
saves a large amount of processing power. We implemented a 
criterion, according to which a line is ignored if its line-core 
opacity is at most 10% larger than the continuum opacity at 
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all depth-points in the atmosphere. The line-core opacity is 
calculated as: 



770 = 4.9947 x 10 



■21 



lQiogte/)( W; _ nu ) 
PA4 



(3) 



where gi is the statistical weight of the lower level, / is the 
oscillator strength, «/ and n u are the atomic lower and upper 
level number densities, p is the gas density and AAd is the 
Doppler width (see Eq @ below). If all these parameters are 
given in c.g.s units, the resulting t/q is in units of cm 2 g -1 . 
Line profiles. Computing Voigt functions (? ?) for a large 
number of lines, wavelength points and many Zeeman com- 
ponents is the dominant fraction of the spectral synthe- 
sis time, even when using an optimized vector routine. 
Extending the Voigt profiles too far into the wings results 
in a lot of unnecessary processing. On the other hand, falling 
short produces inaccuracies near the continuum. A good es- 
timate of the width of the line is therefore very important 
in optimizing the efficiency of the code. We have imple- 
mented an adaptive scheme that starts with a standard width 
of 0.01 nm and extends the interval until the opacity at the 
boundaries is smaller than 1% of the continuum opacity. 
Figure [2] shows the Call line at 854.2 nm in the solar-like 
model. This is perhaps the line that shows the most conspic- 
uous glitch (marked by arrows in the figure) in the transi- 
tion from line to continuum. The glitch is noticeable only 
marginally in Stokes / and is negligible in Stokes V. 
Wavelength sampling. Since the Doppler width of the lines 
increases with wavelength (and the Zeeman splitting in- 
creases with the square of the wavelength), one needs more 
resolution in the blue than in the red side of the spectrum. 
We divided the spectra in 100 nm intervals and used a sam- 
pling typical of solar observations in the visible (^20 mA at 
650 nm). This sampling is then scaled linearly with wave- 
length. It may be seen that even for cool atmospheres and 
relatively heavy ions our sampling is still smaller than the 
Doppler width of the lines, given by: 



AA r 



2kT v mk 
+ 

niA c 



(4) 



where k is the Boltzmann constant, T is the temperature in 
the atmosphere, is the atomic mass, v ml£ . is the microtur- 
bulence and c is the speed of light. 

Weakly Polarizing Media. In many cases, the radiation field 
is weakly linearly polarized, while r\Qjjlr\i <S 1, and 
magneto-optical terms are negligible (pQjj.vMi !)■ I n 
that case, the coupled system of Eq (Q]) simply reduces to 
two decoupled differential equations for I ±V that are easily 
integrated. If circular polarization is also weak, we can ne- 
glect the -J]v V from the equation for the intensity (this is the 
Weakly Polarizing Media approximation of ? ?) and we get 
an independent equation for the intensity /, and an easily in- 
tegrable equation for V. If the criterion above is satisfied for 
all depth-points in the atmosphere, the code automatically 
switches to this method for the formal solution. 



3. Discussion 

Because the Zeeman splitting increases with wavelength, one 
might expect infrared regions to be more promising targets for 
magnetic field observations. However, this is not necessarily the 
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Fig. 1. Blue continuum spectrum for a solar-like model. The 
lower panel shows a close-up of the region marked by vertical 
dotted lines in the upper panel, rescaled to show the discretiza- 
tion of the continuum. Notice, however, in the ordinates scale 
that this effect is very small. 
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Fig. 2. Spectral profile of the Ca II line in the solar-like model, 
showing some small glitches (marked by the arrows) in the 
line-to-continuum transition of Stokes / (upper panel). Stokes V 
(lower panel) is much less sensitive to this effect. 



case. Figure [3] shows the total amount of polarization (in abso- 
lute value) measured in each 100 nm interval of the solar- like 

model spectrum, i.e., f 1 \V(A)\/I c (A)dA. By far, the blue side of 
the spectrum carries a much higher density of polarization sig- 
nal than the red or infrared. In fact, for the strongest fields, the 
region from 400 to 500 nm exhibits nearly as much signal as the 
rest of the spectrum combined (from 500 to 1600 nm). The rea- 
son for this is two-fold. First, the blue region is heavily crowded 
with a large number of closely packed spectral lines. This is ev- 
ident from Table [T] where we can see that the number of lines 
per wavelength interval decreases monotonically towards longer 
wavelengths. Second, blue lines are typically much stronger than 
red or infrared lines. The degree of polarization produced by a 
Zeeman-split line depends not only on the magnetic field and the 
component splitting, but also on the line strength. 

We can see in Figure [3] the balance among all the different 
effects for this model. The amount of Stokes V signal is largest 
at 400 nm and decreases monotonically across the visible spec- 
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Fig. 3. Integrated absolute value of Stokes V (in units of the 
continuum intensity) over 100 nm bins of the solar-like model. 
Solid: 100 G. Dashed: 1000 G. Dotted: 5000 G. The scale of the 
dashed (dotted) curve has been reduced by a factor of 10 (30) in 
order for all three curves to be plotted together. 



trum. A minimum is reached between 1200 and 1300 nm, after 
which the three curves start to increase again all the way up to 
1600 nm. For the stronger fields, the signal drop with wavelength 
is considerably more pronounced and the increase after 1300 nm 
is much weaker. This is because the amount of polarization pro- 
duced by a line increases with the field strength only in the so- 
called weak-field regime (see, e.g., ? ?). This is the regime in 
which the Zeeman splitting is smaller than the Doppler width of 
the line. If all other parameters are constant, the ratio of Zeeman 
splitting to Doppler width increases linearly with wavelength. 
As we increase the field strength, red and infrared lines are more 
likely to depart from the weak-field regime before blue lines do. 
In the saturation regime, the polarization amplitude does not in- 
crease with the magnetic field any more. In summary, observing 
longer wavelengths is only advantageous (in terms of polariza- 
tion amplitude) for relatively weak fields (~100 G). 

If the magnetic field is not homogeneous and covers only 
a fraction of the stellar surface, then the Stokes V spectrum 
will have its overall amplitude decreased according to the filling 
factor (fractional surface area occupied by the magnetic field). 
Generally, it is not straightforward to separate the field strength 
from its filling factor. Figure[3]gives us a clue on how this might 
be accomplished, since the shapes of the curves for different in- 
trinsic field strengths are clearly different. If one were to take the 
ratio of polarization signal in the 800 nm to the 400 nm band, this 
ratio would be larger for the weaker fields and decrease as the 
field becomes stronger. This is a similar concept to the line-ratio 
technique used by ? (?) to propose that the solar magnetic net- 
work is actually composed of strong kG fields occupying small 
filling factors. 

There are many other lines that could potentially be used to 
disentangle intrinsic field strength from filling factor. Table [2] 
lists a number of spectral lines that have the same line opac- 
ity (same lower level excitation potential and log(g/)). The Fe I 
lines at 412.180 and 899.956 nm are particularly promising (see 
? ?), although the disparity in wavelength may introduce com- 
plications (most notably, the background opacities are sensibly 
different). There is a number of Fel, Cal, Til and Si I that also 
have the same Lande factor, and are therefore virtually identical 
except for their wavelength. These lines may be useful, for ex- 




Fig. 4. Stokes V signal in the region around 449 nm. This is 
where the strongest polarization is found for the 5000 G field. 
The vertical dotted line marks the central wavelength of a strong 
Fe I transition at 448.9739 nm. The half-amplitude of this profile 
is approximately 44% of the continuum intensity. 

ample to determine accurate atmospheric temperatures through 
the calculation of their equivalent width ratio. 

Table[3]lists the wavelengths with strongest polarization sig- 
nal in the solar-like model. We can see that for weak fields the 
strongest signals are in the infrared. As the field strength in- 
creases, most of the infrared lines saturate and the strongest 
signals are observed in shorter wavelengths. Note that it is not 
straightforward in general to associate specific spectral lines to 
those wavelengths because some times (especially in the blue), 
the polarization signal arises from a blend of several lines. For 
example, the strongest signal in the solar spectrum produced by 
weak fields (~100 G) comes from a Fel line with central wave- 
length at 1 159.358 nm with an effective Lande factor of 2.5. For 
very strong fields, on the other hand, the strongest signal comes 
from a very crowded region around 448.97 nm. Figure [4] shows 
the Stokes V profile in this region. The most conspicuous feature 
is produced mostly by an Fe I line at 448.9739 nm (vertical dot- 
ted line), but there are other lines with very close wavelengths 
that contribute to it. 

The variation with spectral type is shown in Fig [5] for 
log(g)=4.5. The later spectral types (left of the figure) have a 
very prominent blue side of the spectrum, with a steep decrease 
towards the red. Earlier types exhibit a rather flat spectral depen- 
dence, with a minimum around 700 nm. Notice also that there 
is an overall decrease of the polarization signal for hotter stars. 
This is due to a combination of three different physical effects. 
First, the thermal broadening of spectral lines increases with 
T e ff, smearing out the polarization signal. Second, the increase 
in excitation and ionization state of the atoms in the plasma re- 
duces the number and strength of spectral lines (especially, the 
abundant low-excitation Fe 1 lines). This effect is particularly 
important in the blue. Finally, if we compute the source function 
(S ) and its gradient {dS /dr) and plot the ratio (dS / c/r) /S in the 
photosphere as a function of T e jf, we find that this ratio has a 
maximum around 5000 K and then decreases monotonically for 
hotter stars. 

4. Conclusions 

Stellar polarimetry is becoming increasingly important in 
Astrophysics and it is now possible to obtain polarized spectra 
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Fig. 5. Integrated absolute value of Stokes V (in units of the continuum intensity) over 100 nm bins for varying spectral types and a 
lkG field. 



over a broad range of wavelengths (examples are the Narval and 
Musicos instruments at Pic du Midi, ? ?; or ESPaDOnS at the 
Canada-France-Hawaii Telescope, ? ?), or over a narrower range 
with very high spatial resolution and large aperture (FORS at the 
VLT; e.g., ??;??). 

We have used a Stokes synthesis code and a large array of 
LTE models to produce an atlas of the polarization signal ex- 
pected from a wide variety of magnetic stars. We hope that this 
tool will be useful for the planning of future polarimetric observ- 
ing runs and perhaps for the interpretation of the data, as well. 

For example, it may be possible to (coarsely) generalize the 
popular least-squares deconvolution technique (? ?) to inter- 
mediate and strong magnetic regimes. In order to do this, one 
would take the atlas for the specific star under consideration and 
add all the spectral lines together, as is usually done when ap- 
plying least-squares deconvolution, for different magnetic field 
strengths. By comparing the observed profile to those obtained 
from the atlas, we may be able to estimate the stellar magnetic 
field. 

In principle, it is conceivable to use the synthesis module de- 
veloped here for the application of least-squares inversions (e.g., 
??;??) to stellar spectro-polarimetric observations. It would 
be a siginifcant computational undertaking, however, since this 
kind of inversions typically require of the order of ~100 syn- 
theses (or, if one wishes to use numerical response functions, 
this figure would have to be multiplied by the number of free 



parameters). An interesting alternative may be the use of inver- 
sions based on Principal Component Analysis combined with a 
look-up algorithm (? ?;??;? ?). To this aim, it would be desir- 
able to develop a more refined grid with a better sampling of the 
magnetic field and, possibly, a variety of geometries. 

A final caveat in using our atlas is that it has been derived 
from a very simple magnetic toplogy. Active stars are likely to 
exhibit more complex magnetic structuring. The integration of 
different field orientations over the stellar disk would result in an 
overall decrease of the circular polarization signal with respect 
to the synthetic profiles presented here. However, the relative 
scaling between different lines and spectral regions should be 
approximately correct. 

Acknowledgements. This work has been partially funded by the Spanish 
Ministerio de Educacion y Ciencia through project AYA2004-05792 
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Table 2. Compilation of spectral line pairs with identical line 
opacities 
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Wavelength 


Excitation 
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(nm) 
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Fel 


400.560 


4.186 


-2.131 


2.0000 


Fel 


453.894 


4.186 


-2.131 


0.8333 


Fel 


402.473 


3.241 


-0.710 


0.6250 


Fel 


408.530 


3.241 


-0.710 


1.7083 


Til 


406.509 


1.053 


-0.760 


1.5000 


Til 


469.876 


1.053 


-0.760 


1.0000 


VI 


410.516 


0.267 


-0.130 


1.5000 


VI 


440.851 


0.267 


-0.130 


1.4667 


Fel 


412.180 


2.832 


-1.300 


0.6667 


Fel 


899.956 


2.832 


-1.300 


1.5000 


Nal 


421.606 


2.104 


-2.350 


1.1000 


Nal 


454.518 


2.104 


-2.350 


1.1667 


Fel 


428.582 


3.640 


-1.780 


1.1667 


Fel 


496.870 


3.640 


-1.780 


0.8333 


Cal 


435.508 


2.709 


-0.430 


1.0000 


Cal 


452.693 


2.709 


-0.430 


1.0000 



Table 3. Wavelengths (in nm) exhibiting strongest polarization 
in the solar-like model. The numbers in parenthesis give the de- 
gree of polarization (relative to the local continuum). 



100 G 1000 G 5000 G 



1159.34 (7.18%) 


525. 


.02 


(35 


.38%) 


448. 


.97 


(44. 


.36%) 


846.83 (6.56%) 


423. 


.27 


(34 


.63%) 


420. 


.66 


(44 


.13%) 


1449.77 (6.27%) 


425. 


.83 


(34 


.21%) 


411. 


.17 


(43. 


.63%) 


525.02 (5.94%) 


524. 


.75 


(33 


.64%) 


431. 


.28 


(43. 


.37%) 


1564.84 (5.78%) 


408. 


.08 


(33 


.62%) 


400. 


.16 


(43. 


.24%) 


1386.41 (5.74%) 


513 


.14 


(33 


.36%) 


417. 


.49 


(43. 


.11%) 


1168.98 (5.71%) 


461. 


.33 


(33 


.33%) 


408. 


.29 


(43. 


.02%) 


1142.22 (5.45%) 


454. 


.59 


(31 


.96%) 


412. 


.12 


(42. 


.67%) 


1456.57 (5.45%) 


550 


.67 


(31 


.95%) 


429. 


.14 


(42. 


.65%) 


513.14 (5.38%) 


522. 


.55 


(31 


.89%) 


457. 


10 


(42. 


.55%) 


524.75 (5.37%) 


448. 
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(31 


.66%) 
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44 


(42. 


.55%) 


525.01 (5.26%) 


429. 


.14 


(31 


.10%) 
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.92 


(42. 


.53%) 


1399.74 (5.25%) 


465. 


.44 


(30 


.81%) 


405. 


.55 


(42. 


.43%) 



Fel 443.834 3.686 ~ -1.630 2.5000 

Fel 466.806 3.686 -1.630 1.7500 

Mnl 445.301 2.941 -0.490 1.5000 

Mnl 473.909 2.941 -0.490 0.8000 

Fel 456.531 3.274 -2.510 2.0000 

Fel 479.125 3.274 -2.510 0.5000 

Fel 457.982 3.071 -2.830 1.7500 

Fel 718.915 3.071 -2.830 1.0000 

Fel 461.929 1603 Tl20 1.8333 

Fel 463.801 3.603 -1.120 1.7917 

Fel 463.585 2^845 ^420 2.2500 

Fel 635.503 2.845 -2.420 1.0000 

Til 469.133 1.067 -0.840 1.1667 

Til 586.645 1.067 -0.840 1.1667 

Fel 480.065 4.143 -1.260 1.0000 

Fel 518.792 4.143 -1.260 1.0000 

Sil 540.466 5.614 -0.530 2.0000 

Sil 548.898 5.614 -0.530 2.0000 

Sil 540.466 5.614 -0.530 2.0000 

Sil 623.732 5.614 -0.530 2.0000 

Sil 540.985 5.616 -0.690 2.0000 

Sil 624.447 5.616 -0.690 2.0000 

Sil 548.898 5.614 -0.530 2.0000 

Sil 623.732 5.614 -0.530 2.0000 

Fel 560.294 1430 ^0851 0.7500 

Fel 565.853 3.430 -0.851 2.2500 



Sil 565.492 5.614 -0.930 2.0000 

Sil 612.502 5.614 -0.930 2.0000 



